#include <spu_mfcio.h>

//double myrberg

void recu(vector float * x,  vector float * y, vector float a, vector float b)
{
  *x = spu_mul(*x, *x);
  *x = spu_mul(*x, a);
  *x = spu_add(*x, b);
  *y = spu_mul(*y,*y);
  *y = spu_mul(*y, a);
  *y = spu_add(*y, b);	
}
